Fish-tail and vortex dynamics in Ni doped iron-pnictide BaFei,s2Nio.i8As2 



O 

(N 



o 
o 

u 

Oh: 



O 

> 
00 

O 



X 



S. Salem-Sugui Jr./ L. Ghivelder,^ A. D. Alvarenga,^ L.F. Cohen, Huiqian Luo,* and Xingye Lu** 

' Instituto de Fisica, Universidade Federal do Rio de Janeiro, 21941-972 Rio de Janeiro, RJ, Brazil 
"^Instituto Nacional de Metrologia Normalizagdo e Qualidade Industrial, 25250-020 Duque de Caxias, RJ, Brazil. 
"'The Blackett Laboratory, Physics Department, Imperial College London, London SW7 2AZ, United Kingdom 
'^National Lab for Superconductivity, Institute of Physics and National Lab 
for Condensed Matter Physics, P. O. Box 603 Beijing, 100190, P. R. China 

(Dated: June 15, 2011) 

We study the vortex-dynamics of a BaFe\,s2Nio.-if!,As2 crystal with Tc = 8 K, by measuring flux- 
creep over the second magnetisation (or fish-tail) peak for both H|jc-axis and H||ab planes. Magnetic 
relaxation data show an anomalously long initial stage of relaxation, lasting approximatelylO min- 
utes for H||c-axis and 2-3 min for H||ab, resembling a transient effect with a lower relaxation rate, 
which is followed by the usual log(time) relaxation. Interestingly, study of the relaxation rate R 
vs H for both stages of relaxation and for both field directions, are featureless over the field range 
associated with the fish-tail. The same trend was confirmed by plotting R vs T obtained from 
flux-creep data measured as a function of temperature for a fixed field (H||c-axis). A plot of the 
activation energy U(M,T) calculated from the time relaxation of the magnetisation at a fixed field 
also shows a smooth behavior further supporting the view that the fish-tail peak is not associated 
with a crossover in vortex pinning regime within the collective pinning scenario. 

PACS numbers: 74.70.Xa,74.25.Uv,74.25.Wx,74.25.Sv 



The observation and study of the second magnetisation 
or fish-tail peak in pnictidesi has attracted the atten- 
tion of an increasing scientific audience. The relatively 
large flux-creep observed in this material system allows 
detailed studies of the vortex-dynamics and comparison 
to a rich variety of models^.. Untill now the fish-tail effect 
has been studied in SmFeAsO^.i^FQ i^ , NdFeAsOo.sd'^ , 
BaFe2-xKxAs-2r^, BaFe2-xCoxAs^^—, and more re- 
cently in LiFeAiM and PrFeAsOo.eoFOAi^. Different 
mechanisms have been claimed to account for the fish- 
tail effect in these compounds (see for instance a table 
presented in RefH) evidencing the need to build up a 
more comprehensive picture through additional studies, 
particularly on systems not yet investigated. 
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In this work, we present a vortex-dynamics study per- 
formed on the recent synthesized Ni doped iron pnictide 
superconductor system BaFei,s2Nio,isAs 
an overdoped single crystal with x=0.18, 
K), ni=78 mg, and average dimensions 1.1x0.4x0.025 cm, 
for magnetic fields applied along the ab-planes direction. 
In this geometry, a possible misalignment of the sample 
with respect to the direction of the applied magnetic field 
is estimated to be smaller than 2 degrees, which assures 
that any contribution of the c-axis component to the 
magnetization measured with II|| ab-planes is neglegible. 
A small piece with m=23.4 mg and average dimensions 
04x0. 4x0. 02cm was broken from the large sample for the 
measurements with H||c-axis. Details on sample prepa- 
rations and physical properties of the large high-quality 
single crystals can be found in RefJ5.. Isothermal M(II) 
curves for the crystals for H|| ab-planes and H|| c-axis ex- 
hibit the fish-tail peak, the maximum occurring at an 
applied field Hp. The vortex dynamic investigation was 
performed by measuring isofield and isothermic magnetic 
relaxation, M vs time curves, over a period of 1 hour for 



magnetic field values lying below and above Hp on both 
branches of the M(H) curves. All data was obtained after 
cooling the system from above Tc to the desired temper- 
ature in zero applied magnetic field (but in the presence 
of the earth magnetic field), which is called the zfc proce- 
dure. Magnetization data was obtained using commercial 
magnetometers based on a superconducting quantum in- 
terference device (SQUID) (MPMS-5T Quantum Design, 
and Criogenics-6T) . The charge and discharge rate of the 
magnet were set equal to 100 Oe/s during the experi- 
ment. 

Figure 1 shows isothermal M(H) curves exhibiting the 
fish-tail effect, as obtained for both geometries. M(H) 
curves at 3K and 4K for H||c and at 3K for H||ab are 
plotted with correspondent magnetic relaxation data ob- 
tained on both branches. We refer to this as the isother- 
mal method. The dotted line located at H=4.6 kOe in 
Fig. la show that once this field value is fixed it is pos- 
sible to go from below Hp to above it by increasing the 
temperature from 3 K. This second method, which we 
refer to as the isofield method, provides an independent 
check of the relaxation data taken using the isothermal 
route and allows the creation of an activation energy 
curve with temperature U(T) as in Ref flTlfTsl . As de- 
picted in Figure 1, all M(H) curves are quite symmetric 
relative to the x-axis evidencing that vortex dynamics 
is mostly due to bulk pinning. Also, since the equilib- 
rium magnetization, Meq—{M^ + M^)/2w0, where M+ 
and M~ refer to the magnetization in the upper and 
lower branches of the hysteresis curve respectively, we 
use M instead M — Meq in all curves. The upper inset of 
Fig.l shows a plot of the critical current Jc vs H at 3K 
for both geometries as estimated from the correspondent 
hysteresis curves shown in Fig.l by using the well know 
expression Jc = 20AM {emu/ cm^)/[a{l - a/36)] where 
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FIG. 1: Isothermic M(H) curves, a) H||c-axis, at 3, 3.5,4, 4.5 
and 5 K. The upper inset shows Jc(T=3K)vs H in logarithmic 
scale for both field directions. Solid lines in the region of the 
peak represents a fit to an expression of Ref[l(^. The lower 
inset shows Hp vs T extracted from M(H[) curves for both 
field directions; b) H||ab-planes. 



b>ai2,. The lower inset of Fig.l shows a plot of the fish- 
tail peak field. Hp, as a function of temperature extracted 
from the curves of Fig.l for both field directions, where 
dotted lines are only a guide for the eyes. We observe 
from the inset figures that Jc(Hp||ab)/ Jc (Hp ||c)w3. 3 and 
Hp(H||ab)/Hp(H||c)w3.6. Although neither of these pa- 
rameters is a fundamental measure of the superconduct- 
ing anisotropy (for that we would have to extract the 
upper and lower critical fields), they clearly are consis- 
tent with each other, and the fact that the anisotropy 
expected for this system is of the order of 3^. 



Figure 2 show selected magnetic relaxation data ob- 
tained on both branches of the M(H) curve at 3 K for 
H||c-axis. These curves exemplify the general behavior 
observed on all relaxation curves obtained in this work, 
namely, all curves show first a slow relaxation rate, and 
then an increase in the relaxation rate after a time to~10 
minutes after relaxation starts for H||c and, to«2-3 min- 
utes for H||ab. This anomalous behavior showing two 
distinct linear behavior of M with log(t) or two distinct 
time windows, can not be explained in terms of a pin- 
ning crossover, since in that case the crossover always 
occur from a higher rate of the initial stage of relax- 
ation (corresponding to flux jumping over low activation 
energy pinning sites) to a lower rate (corresponding to 
fiux jumping over higher activation energy pinning sites), 
as occurring in the well known surface to bulk pinning 
crossover^". This scenario is the opposite of what is ob- 
served here. There are two other possible explanations 
for the non-linearity of M(t) vs logt curves. The first is 
due to an eddy current induced on the sample by the 
ramp rate dH/dt while charging the magnet^!. After the 
magnet is charged and dH/dt=0 the induced eddy cur- 
rent decays due to flux-creep producing a transient region 
which possesses a lower magnetic relaxation rate when 
compared to that due to the bulk pinning. The duration 
time of this transient relaxation is inversely proportional 
to the charge magnet rate dH/dt and it is expected to 
?-i^S^xj^§^««itxjj»>t)«i«^ I last less than a minute for a rate dH/dtwlOO Oe/s^^ as 
ii^-*? used here. This effect referred to as an initial settle time, 
is commonly observed in M(t) curves, and is inversely 
proportional do the ramp rate of the magnetic fleld^i 
and found to be of the order of 1 minute usually. It 
should be mentioned that the ratio between the times 
To(H||c)/To(H||ab)ss4 for our data qualitatively agrees 
with the hypothesis that the anomalous initial decay of 
M(t) is related to the transient relaxation associated to 
dH/dt, once r^KJia/ {dB / dt^ where Ji represents the 
eddy current (J(H||ab)«3J(H||c see Fig. 1), a is the sam- 
ple thickness perpendicular to H (a=0. 025cm for H||ab 
and 0.4cm for H||c), and dB/dt is the rate of the mag- 
netic induction field B. The other possible explanation, is 
based on the existence of twin-boundaries^^ which may 
allow a continuous increase of vortices pinned between 
twin-boundaries (appearing in the first time window). 
As a result, after the density of vortices pinned by twin- 
boundaries increases above a certain value, the pressure 
overcomes the barrier producing an increase in the rate 
of relaxation at the second time windowS^.. This expla- 
nation would apply as well to our experiment, since the 
efficiency of twin-boundaries, as pinning sites, changes 
and become less effective as the sample rotates from H||c 
to H||ab (this because many TB aligned parallel to the 
fied when H||c can become perpendicular to the field after 
the 90 degree rotation). However, after careful observa- 
tion of the sample surfaces in a polarized light micro- 
scope, there are no visible twin-boundaries, and up to 
this point, twinning has only been observed in under- 
doped not overdoped crystals. Consequently we assign 
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FIG. 2: Selected M vs time curves for H||c: a) lower branch of 
M(H). b) upper branch of M(H). Solid lines are only a guide 
to the yes evidencing the increasing relaxation rate for t>ro. 
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FIG. 3: Plots of the relaxation rates; R vs H for both branches 
of M(H) curves for: (a) and (b) H||c-axis; (c) H||ab; (d) R vs 
T for H=4.6 kOe, Hllc. 



the initial relaxation to a transient effect. 

We analyze flux-creep data by obtaining the relaxation 
rate R = dM/ dint of the second stage of relaxation (cor- 
responding to the region t>ro) for each M(t) curved"—. 
Figure 3 show the results from relaxation obtained in the 
lower branch (i?mcr) and in the upper branch (Hdecr) 
of M(H) curves for both geometries. We mention that 
we also obtained values of R(t<r) for the first stage of 
relaxation, and observe that plots of R(t<r)vs H ap- 
proximately follow the same trend as the plots in Fig. 3. 
This fact may evidence that the same major bulk pinning 
mechanism is present in both stages of relaxation, sup- 
porting our hypothesis that the anomalous first stage of 
relaxation is related to a transient effect rather than to 
a twin-boundary assisted effect. The main information 
one can extract from the many plots of Fig. 3 is the ab- 
sence of any feature, minimum or maximum, located near 
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FIG. 4: The scaled activation energy U(M,T) for H=4.6 kOe 
are plotted against M. The solid line is a logM fit. 



the correspondent fish-tail peak position represented by 
a vertical arrow in each figure. This fact strongly sug- 
gests that the fish-tail peak observed in M(H) curves for 
both field directions is not due to a pinning crossover. 

We also obtain flux-creep data from the isofield 
method, with a path represented in Fig. la by a verti- 
cal dotted line for H=4.6 kOe. M(t) data were obtained 
by reaching the desired temperature in a zfc procedure, 
followed by an increase of the field up to 4.6 kOe and 
measuring M(t) for 60 min. Results of R=dM/dlnt ob- 
tained this way for all M(t) data is presented in Fig. 
3d. As in other plots of R vs H shown in Fig. 3, the 
plot of Fig. 3d does not show any marked feature as Hp 
(located at T=4K) is crossed. The isofield activation en- 
ergy curve U(M,T), or U{M,T) / g{T /T^) where g{T/Tc) 
is some scaling function, is believed to be a smooth func- 
tion of M, where its behavior with M can provide in- 
formation on the pinning mechanisroiiii^'2^. Here we 
performed an analysis of the activation energy as shown 
previously in Ref ilTlilSl using the expression 

U = -Tln{dM{t)/dt) + CT 

where C is a constant which depends on the hoping dis- 
tance of the vortex, the attempt frequency and the sam- 
ple size. The constant C is adjusted in a manner such 
that all U{M,T)/g{T/Tc) curves plotted against M fall 
on a smooth curve. If all data does not fall on one smooth 
curve with a fixed value of C, it is an indication that the 
data does not fall into one pinning regime. Figure 4 show 
the results of the analysis. The appropriate scaling func- 
tion for our case is g(T/Tc)=(l - T/T^f/'^ also used in 
Reffisl and the constant C=14; similar values of the con- 
stant C have been found for high-Tc cupratesiiii^ and 
pnictides®. The plot of U(T) show an almost perfect log- 
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arithmic with M, indeed suggesting that a single pinning 
regime is operating in the temperature range studied at 
the fixed field of 4.6 kOe. 

As a pinning regime crossover does not appear to be 
an appropriate explanation for the fishtail effect it is 
interesting to explore whether the data can be under- 
stood within the model proposed for a vortex lattice 
phase transition^^, using the expression Jc{B)=A/[{B — 
BpY + (A_B)^]^/^ where A is a fitting parameter and 
AB is the peak width. The M vs H curves have been fit 
to this expression, converting the magnetisation data to 
critical current as described previously and using B«H. 
The good-quality fittings conducted on Jc(T=3K) for 
both field directions are shown as solid lines in the 
curves of the upper inset of Fig. la. By considering 
Jc=AM(emu/cm'^) the values of the fitting parameters 
are: A=l.lxlO'* G^/^ g^^^ AB=U.2 kOe for H||c, and 
A=4.4xlO^ G^/^ and AB=A6 kOe for H||ab. Note that 
the fittings produced larger values for the peak width 
compared to values obtained in ReflT2l. It is intersting 
that the data fits well to this expression, although we 



recognize that this in itself is not conclusive evidence of 
a phase transition in the vortex lattice. Moreover we 
would expect a change in vortex dynamics as a signa- 
ture that a phase transition has occurred, and we see no 
evidence for this. 

In conclusion, our experiments show that the fish- 
tail peak appearing on M(H) curves for H||c and II||ab 
in Bai^ei. 82-^*0. i8^'S2 does not appear to be associated 
with a softening in vortex pinning prior to melting, nor 
a change of pinning regime within a collective pinning 
model scenario. We observe an anomalously long relax- 
ation period at the initial stage of relaxation lasting wlO 
minutes for II||c which is likely to be a transient effect 
due to the field rate dH/dt while charging (discharging) 
the magnet. This transient effect appears to be a feature 
intrinsic to the BaFeNiAs system we have studied, as we 
have also found it present in preliminary data obtained in 
a second crystal with x=0.1 (T'c«20K) for II||ab-planes. 
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